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ABSTRACT 

We made a phase-resolved timing analysis of GRS 1915+105 in its p state and ob- 
tained detailed p cycle evolutions of the frequency, the amplitude and the coherence of 
low-frequency quasi-periodic oscillation (LFQPO). We combined our timing results 
with the spectral study by Neilsen et al. to perform an elaborate comparison analy- 
sis. Our analyses show that the LFQPO frequency does not scale with the inner disk 
radius, but it is related to the spectral index, indicating a possible correlation between 
the LFQPOs and the corona. The LFQPO amplitude spectrum and other results are 
naturally explained by tying the LFQPO to the corona. The similarities of the spec- 
tra of variability parameters between the LFQPO from p state and those from more 
steady states indicate that the LFQPOs of GRS 1915+105 in very different states seem 
to share the same origin. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: individual 
(GRS 1915+105) 
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1. Introduction 



GRS 1915+105 is a binary system that was discovered by WATCH on board GRANAT 



in 1992 (Castro-Tirado et al. 



about 11 kpc (e.g.. 



(|Zhang et al. 



1997 



1992). 



? ender et al. 



1999; 



McClintock et al 



t is located in our Ga laxy at an estimated distance of 



Zdziarski et al. 



with mass 0.8 ± 0.5 M Q as the donor (IHarlaftis & Greiner 



20051) . containing a spinning black hole 



2006) with mass 14 ± 4 M ? , and a K-M III gia nt star 



2004 



Greiner et al. 



2001b). The 



orbital separatio n of the binary components is about 108 ± 4 R Q and the orbital period is 



33.5 ±1.5 days ( Greiner et al. 



superluminal radio jets ([Mirabel & Rodriguez 



2001 a]). As the first microquasar found, GRS 1915+105 produces 



1994 



Fender et al. 



1999). The count rate and color 



characteristics are extremely complex and the light curves of the source are usually classified 
into 12 v ariability classes wh ich are regarded as transitions between three basic states, states A, 



B and C (Belloni et al 



1 



2000). Lots of low-frequency (~0.5-10 Hz) quasi-periodic oscillations 



(LFQPOs) are found in GRS 1915+105, providing us ideal subjects for the study of LFQPO. 
Fruitful results have been obtained from the LFQPO of GRS 1915+105. It is revealed that 



the sr. 


>ectra, i.e., the thermal and power-law components or the source intensity 


(e.g.. 


Chen et al. 


1997: 


Markwardt et al. 


1999: Muno et al. 


1999.2001: 


Trudolvubov et al. 1999: 


Reig et al. 


2000: 


Tomsick & Kaaret 


2001 


). The LFQPO amplitude is inversely correlated with the source flux or 



LFQPO frequency (e.g. 



Muno et al. 


1999; 


Reig et al. 


2000; 



Trudolvubov et al. 



1999). As the 



LFQPO frequency increases, the tem perature of the inner accretion disk increase s and the radius 



of the inner accretion disk decreases (Muno et al 



1999 



Rodriguez et al. 



2002b|). The LFQPO 



amplitude i ncreases with photon energy and in some cases i t turns over in the higher en ergy 



bands (e.g., 



Tomsick & Kaaret 



2001 



Rodriguez et al. 



2002a, 



2004; 



Zdziarski et al 



2005). As 



the centroid frequency of the LFQPO increases, the relation betwee n the LFQPO fr equency and 



photon energy evolves from a negative correlation to a positive one (lOu et al 



2010). In addition, 
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three additional combined patterns of the negative correlation and the positive one were disco vered 



in a systematic study of the LFQPO frequency— photon energy relationship (|Yan et al 



2012 ) 



Although the results mentioned above enable a good understanding of LFQPO phenomenol- 
ogy, we are puzzled by the ambiguity that the LFQPO is correlated with b oth the corona/jet and 



Neilsen et al. 



(1201 lL hereafter NRL11) 



accretion disk. The LFQPO origin is still a mystery, 
investigated the physical changes and the LFQPO evolution of the p variability in GRS 1915+105 
through a phase-resolved spectral and timing analysis. In order to reveal more clues about the 
origin of the LFQPO and more details about evolution of the p cycle, in this work we make 
a detailed comparison between the results of spectral and timing analyses and investigate the 
linkage between the p-associated LFQPOs and the LFQPOs from GRS 1915+105 during more 
steady conditions. In addition, we attempt to use the ratio of the power density spectrum (PDS) 
continuum amplitude to the LFQPO amplitude to track the origin of aperiodic X-ray variability 
from accretion flows. 

The observation and data reduction methods are described in Section 2, the results are 
presented in Section 3, and the discussion and conclusion are given in Section 4. 



2. Observation and Data Reduction 

With a method simialr to that of NRL1 1 , we produce the average phase-folded light curve for 
60405-01-02-00, the same RXTE observation of GRS 1915+105 that was analyzed by NRL11. 

We perform our analyses with the HEASOFT version 6.10 package. The light curves 
with a time resolution of 1 s are extracted from the binned mode data (B_8ms_16A_0_35_H_4P) 
in 2.0-14.8 keV and the event mode data (E_16us_16B_36_ls) in 14.8-60 keV. The data are 
accumulated from four PCUs in combinations of {0, 2, 3, 4} and {1, 2, 3, 4}. Good time intervals 
are defined as follows: a satellite elevation over the Earth limb > 10° and an offset pointing 
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< 0.02°. The rows with non-common time are deleted from the two light curve FITS files. Then, 
we add the two dead-time-corrected and background-subtracted light curves and obtain a light 
curve that will be further barycentered. 

We obtain a series of cycles from the barycentered light curve and then cross -correlate them 
with a representative cycle. Maxima in the normalized cross -correlation values correspond to the 
peaks in the light curve. We smooth the cross-correlation with a Gaussian of FWHM 3 s and 
measure the time of maxima via parabolic interpolation. The resulting maxima are phase = 
times which are used to create the average phase-folded light curve. In order to mitigate the effect 
of the initial template choice, we use the folded light curve as a new template and cross-correlate 
it with the entire set of cycles and obtain the final set of = times. 

For the timing analysis, the light curve in a certain phase interval is extracted with a time 
resolution of 8 ms. The PDS is produced wit h the dead-time-corrected Pois son noise level 



subtra cted from the PDS (Mo rgan et al. 



Miyamoto et al. 



1997) and with the normalization of 
(|1992|) . which gives the periodogram in units of (rms/mean) 2 /Hz. The PDS for each 0.02 phase 
and 0.4-0.74 phase are computed on 2 s and 4 s sampling durations, respectively. Following 
Belloni et al.1 (120021) . the PDS is fitted with a model that includes several Lorentzians to represent 



the continuum, the QPOs, and other broad features, respectively. The fitting is limited over the 
0.5-20 Hz frequency band. 

In order to study the spectra of the variability parameters (e.g., continuum amplitude, LFQPO 
amplitude, LFQPO frequency), we produce a PDS for several energy bands (see Figure |3]). 
The continuum amplitudes in several frequency bands are obtained by subtracting the LFQPO 
contributions (the LFQPO parameters are obtained from t he fitting mentioned above) from the 
full integrals in corresponding frequency bands (see, e.g. JVaughan et al.ll2003|) . Due to the high 
energy dependence of the background, when studying the energy-dependent PDS, we correct the 
continuum amplitude and LFQPO amplitude for background following A nct = A I3iW x with 
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A being the amplitude, S the sourc e net rate, and B the background rate (|Berger & van der Klis 



1994 



Rodriguez & Varniere 



201 lb . The errors are derived by varying the parameters until 



Ax 2 — 1, at 1 a level. 



3. Results 

3.1. Timing Analysis for 0.02 Phase Intervals 

The PDS for each 0.02 phase interval is computed in the 2.0-37.8 keV band and the results 
are shown in Figure [TJ The phase-folded PCA p class light curve (Figure [Eta)) and the shape of 
the PDS power integrated over 0.5-10 Hz (hereafter referred to as the 0.5-10 Hz amplitude, and 
similarly for other amplitudes) (Figure [TJb)) are similar to those of NRL1 1, except that our count 
rate is higher and the 0.5-10 Hz amplitude (2 s sample duration) is a bit lower than NRL1 l's 
2-10 Hz amplitude (1 s sample duration). Based on the behavior of the LFQPO amplitude (Figure 
□Id)), the cycle phase is divided into six intervals (I: 0.02-0.12, II: 0.12-0.26, III: 0.26-0.4, IV: 
0.4-0.74, V: 0.74-0.92, and VI: 0.92-0.02). In interval I, there is no obvious LFQPO. In intervals 
II and VI, the LFQPO amplitude is positively correlated with the LFQPO frequency. Meanhile for 
intervals III, IV and V, the LFQPO amplitude is negatively correlated with the LFQPO frequency. 
In view of the overall situation, the LFQPO amplitude decreases very slightly in interval IV but 
dips in phase intervals 0.12-0.4 (II and III) and 0.74-0.02 (V and VI). The LFQPO amplitude 
in the phase interval 0.12-0.14 is relatively small. As phase increases, the LFQPO frequency 
decreases (II and III), flattens (IV) and then increases again (V). It is very interesting that the 
LFQPO frequency remains constant in interval IV, and drops at <\> = 0.92. We double the LFQPO 
frequencies in interval VI and find that they (the gray points) and the other LFQPO points together 
show a much more symmetric time evolution over the p cycle (Figure HJc)). As phase increases, 
the coherence of LFQPO reveals a rapid increase (II), and then a quick decrease (III) followed by 
a continuous gentle decrease (IV, V, and VI) (Figure (He)). 
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The LFQPO absolute a mplitude is estimated by multiplying th e 



source count rate (see, e.g,. IMendez et al 



19971: 



Gilfanov et al. 



2003 



LFQ PO amplitude a nd 



Zdziarski et al. 



2005). The 



LFQPO absolute amplitude has a behavior similar to that of the LFQPO amplitude. It increases 
very slightly in interval IV but dips in phase intervals 0.12-0.4 (II and III) and 0.74-0.02 (V and 
VI) (Figure [2]). In interval VI, it rises up rapidly. 



3.2. Spectra of Variability Parameters 

Spectra of variability parameters provide us a way to find the correlations between variability 
and the regions where photons with different energy are produced. They are also good indicators 
of where to find the connections between the LFQPOs in the p state and other states. We 
produce the spectra of the 0.5-64 Hz, 3-64 Hz, and 10-64 Hz continuum amplitudes, the LFQPO 
amplitude and the LFQPO frequency for interval IV during which the LFQPO frequency and the 
power-law index are both relatively steady (Figure [Q Figure 7 in NRL1 1). 

The LFQPO frequency almost increases linearly with photon energy when the energy < 15 
keV, and then it levels off (Figure Ha)). The points are fitted with the least squares (see, e.g., 



Greene 



2002). The slope of the best-fit line (red) is 0.0041 ± 0.0010 Hz keV~ x and the adjusted 



R-square is 0.65. 

Figure Sb) shows the 0.5-64 Hz, 3-64 Hz, 10-64 Hz continuum amplitude spectra and the 
LFQPO amplitude spectrum. As the photon energy increases, all these amplitudes increase and 
then flatten out. The ratios of the continuum amplitudes to the LFQPO amplitude are shown in 
Figure He). As photon energy increases, all of the ratios decrease rapidly and then smoothly level 
off. The points can be well fitted with the functio n R(E) = A- DB ln{exp[(£ tr - E) / D] + 1} 



(function (1) in IShaposhnikov & Titarchuk 



2007). We fix the D at 1 keV. The best-fit values for 



the 0.5 - 64 Hz continuum /LFQPO amplitude ratio are A = 1.51 ± 0.08, B = -0.33 ± 0.08 
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keV -1 , and E tT = 7.53 ± 0.95 keV. The fit for the 3-64 Hz continuum/LFQPO amplitude ratio 
gives: A = 1.26 ± 0.09, B = -0.26 ±0.11 keV" 1 , and E tr = 6.99 ± 1.43 keV. The fit for the 
10-64 Hz continuum/LFQPO amplitude ratio gives: A = 1.05 ± 0.10, B = -0.25 ± 0.20 keV -1 , 
and E tr = 6.05 ± 2.20 keV. The errors for the best-fit parameters are standard deviations. 

3.3. Spectral Index-LFQPO Frequency Relation 

In order to further investigate the origin of the LFQPO and the evolution of the p cycle, the 
ratio of the spectral index to the LFQPO frequency is presented in Figure |4] The spectral index is 
obtained from Figure 7 in NRL11. In interval IV, both the LFQPO frequency and the power-law 
index remain almost constant, and the ratio is ~ 0.4. As the LFQPO frequency increases, the 
index increases, while the ratio decreases faster. We must note that the ratio is an artificial quantity 
with arbitrary units and is merely a vehicle to track the correlation between the Compton spectrum 
and the LFQPO frequency versus phase. 

4. Discussion and Conclusion 

In this section, we combine our timing results described in Section 3 with NRLll's spectral 
analysis to investigate the origin of the LFQPO from GRS 1915+105 during its p state. 

In the p state, NRL1 1 showed that the energy spectrum has at least two components: the disk 
emission and the corona emission. The LFQPO may originate from the accretion disk, the corona, 
both the disk and the corona, or the interaction of the disk with the corona. 

As phase increases, the LFQPO frequency decreases (II and III), flattens (IV) and then 
increases again (V) (Figure EDc)), while the accretion disk radius increases continuously over the 
phase range 0.1-0.9, which covers intervals II, III, IV, and V (Figure 7 in NRL1 1). Assuming that 
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the spectral model in NRL1 1 is correct, it turns out that the LFQPO frequency does not scale with 
the inner disk radius, suggesting that the LFQPO cannot be tied to dynamical frequencies in the 
disk. Although the possibility of the LFQPO originating from disk should not be ruled out based 
only on this, it is still a puzzle whether the LFQPO originates from the corona. 

In interval IV which is covered by the slow rise in count rate in NRL11, both the LFQPO 
frequency and power-law index are relatively steady while the radius of the inner disk increases 
significantly (Figure [fl Figure 7 in NRL11). NRL11 argued that the local Eddington limit in 
the inner disk is responsible for the slow rise during which the radius of the inner disk increases 
with luminosity at the constant temperature of the inner disk. The fact that the LFQPO frequency 
and the power-law index remain relatively stable coincidentally indicates a possible correlation 
between the LFQPO and the corona. If the LFQPO originates from the corona, then since the 
LFQPO absolute amplitude in interval IV increases slightly (Figure |2]), we can interpret the slight 
decline in the LFQPO amplitude in interval IV (Figure Ed)) as the result of the increase of the 
disk flux (Figure 9 in NRL11). 

In order to further test whether the LFQPO is from the corona, we investigate the spectra 
of the variability parameters in interval IV. It is found that the LFQPO frequency is positively 
correlated with photon energy (Figure |3ja)) , whi ch is similar to the behavior of some LFQPOs 



from GRS 1915+105 during more steady states (lOu et al. 



2010; 



Yan et al 



2012). It was assumed 



that the inner regions of the accretion flow have a harder 



softer spectrum (e.g., 



Revnivtsev et al. 



1999 : 



spectrum while the outer regi ons have a 



Kotov et al 



2001 



Ingram & Done 



2011 



2012 ). If 



this assumption is correct, then any satisfying model for the LFQPO should predict that the higher 
centroid frequency part of the LFQPO is related to the inner part of the flow, while the lower 
centroid frequency part of the LFQPO is related to the o uter part of flow. The Lense-Thirring 



precession of the hot coronal flow within a truncated disc ([Done et al.1 



2007 



and reference therein) 



can produce the LFQPO, with a frequency that is inversely correlated with the character radius 
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of the hot flow (|lngram et al.ll2009Q . Then, the Lense-Thirring precession model of the LFQPO 
seems to have the potential to explain the observed spectrum of the LFQPO frequency. However, 
although the model could possibly explain the positive correlation between the LFQPO frequency 
and the photon ener gy, it is unclear whether it al one can explain the observed negative and 



intermediate ones (|Ou et al. 



2010 : 



Yan et al. 



2012 



It is also found that as photon energy increases, the LFQPO amplitude increases and then 
flattens (Figure |3£b)). If the LFQPO is related to the corona, then the upward trend in the LFQPO 
amplitude might be a representation of the fractional spectrum of the corona versus the total. The 
flattening of the amplitude spectrum would then reflect the energy where the disk contribution is 
negligible. 

It seems impossible that the LFQPO originates from both the disk and the corona 
simultaneously, due to their essential differences in property and location. Besides, if the LFQPO 
is produced by the interaction of the disk with the corona, then it is hard to explain why the 
frequency of the LFQPO produced by the interaction would increase when the radius of the inner 
disk increases and the frequency of Keplerian motion in the disk decreases. It is therefore very 
likely that the LFQPO is produced in the corona. 



The origin of aperiodic X-ray variability from accretion flows is not 



produced by variations 



Arevalo & Uttley 



Uttley et all 



propagating in the accretion flow (e.g., 



2006 



Titarchuk et al. 



20071: 



Gierlihski et alj 



,yubarskii 



2008 



very 



1997 



clear. It might be 



Kotov et al 



Wilkinson & Uttlev 



2001 



200' 



* 



201 lh . Due to the dissipation or filter effect of the flow, all of the higher frequency 



(e.g., > several Hz) variations created in the outer part of the flow (the disk) may be smoothed out, 



flow (see, e.g., 


ievnivtsev et al. 


1999: 


Nowak et al. 


1999: 


Done e 


tal. 


2007: 


Titarchuk e 


tal. 


2007 


Gierlihski et al. 


2008; 


Wilkinson & Uttlev 


2009[ 


Ingram & Done 


201 


liLeU etalJ 


2011 


). In this 



work, we find that the ratios of the continuum amplitudes to the LFQPO amplitude are roughly 
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constants at photon energy > lOkeV (Figure H c ))> indicating that the continuum in the higher 
energy ban d and the LFQPO m ight originate from the same region of the accretion flow, e.g., 



the corona. 



Ingram et al. 



(|2009|) showed that a hot flow can produce a continuum spectrum and 
a LFQPO simultaneously. As photon energy decreases, the continuum/LFQPO amplitude ratios 
increase rapidly (Figure He)). Assuming that both the continuum in a certain frequency band 
and the LFQPO are produced in the corona and they have the same amplitude ratio in both the 
lower and higher energy bands, this hints that the excess continuum power in the lower energy 
band might be from another part of the flow, e.g., the disk. Then it seems that the disk has 
high-frequency (> 10 Hz) variations (see Figure |3fc) ), which could be explaine d by X-ray heating 



of the disk by the varying corona emission (see, e.g., 



Wilkinson & Uttley 



2009). 



There is no obvious LFQPO in interval I, during which there is a hard pulse identified by 
NRL1 1 . NPvLl 1 argued that some material collides with the hot corona after it has been ejected 
from the inner disk due to disk instability during the hard pulse phase. Based on the presumption 
that the LFQPO is produced in the corona, the absence of the LFQPO may be caused by the 
violent disturbance in the corona owing to the collision. 

The LFQPO frequency decreases smoothly during intervals II and III, which cover NRL1 l's 
hard X-ray tail phase during which a short-lived jet is said to be produced. Then, the jet seems 
to be independent of the LFQPO considering that it has no influence on the LFQPO frequency 
based on the smooth evolution of the LFQPO frequency (Figure [TJc)). However, it seems to 
have reduced the LFQPO amplitude considering the dips in the profiles of the LFQPO amplitude 
(Figure (2d)) and of the LFQPO absolute amplitude (Figure |2]). The decrease of the LFQPO 
amplitude may be caused by the increase in flux of a component which is independent of the 
LFQPO. Besides, it is also possible that the LFQPO amplitude itself has been decreased by a 
certain process, e.g., more accretion material forms the jet but not the corona. 

The LFQPO frequency drops at <fi = 0.92 (Figure [TJc)). NRL1 1 argued that the disk becomes 
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unstable and the inner disk radius decreases rapidly after phase 0.9 due to radiation pressure. 
Apart from the steep rise in the power-law index (Figure 7 in NRL1 1), the drop in the LFQPO 
frequency might be another signal indicating that the corona has been changed significantly. 
However, it is unclear why the LFQPO frequency falls to just half of the hypothetical value. The 
decreases of the LFQPO amplitude and the LFQPO absolute amplitude in interval V might be 
caused by the softening of the source during the corresponding interval (Figure 2e in NRL1 1). 
The increase of the LFQPO amplitude and the steep rise of LFQPO absolute amplitude in interval 
VI might show the impact of disk instability and the quick decrease of the inner disk radius on the 
corona. 

The correlation between the spectral index and the LFQPO frequency, which provides 
us with information ab out the origin of the LFQPO and the evolution of the source, h as 



been well studied (e.g., 



Titarchuk et al 



Shaposhnikov & Titarchu 



2006 



1998 



Vignarca et al. 



2003 



Titarchuk & Fiorito 



2004 



2007|) . It is found that there is a saturation of the spectral 
index for high values of the LFQPO in some black hole binaries. A transition layer model 
was established to explain the LFQPO origin a nd the saturation of the spectral index (e.g., 



Titarchuk et al. 



1998; 



Titarchuk & Fiorito 



2004). For GRS 1915+105 during the y state, the index 



satura tes as the LFQPO frequency is larger than ~ 3 Hz (Figure 1 in 



Shap oshnikov & Titarchuk 



(120071) ). Nevertheless, in the p state, the LFQPO frequency is larger than 6 Hz, while the index 
is positively correla t ed wit h the frequency, although there is a trend of saturation (Figure H). 
Titarchuk & Fiorito! (|2004|) explained a similar phenomenon and argued that it occurs when 
a relatively cold outflow from winds downscatters corona photons and softens the spectrum. 
Considering that it is hard to tell whether the change in spectral index is caused by the corona 
change or the wind change, we should be cautious about adopting this argument to interpret 
our result, though there are some agreements between the evolution of the flux of the Fe XXVI 
absorbed line and the evolution of the ratio of the spectral index to the LFQPO frequency (Figure 
12 in NRL11; Figure gj. 
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In interval IV, the LFQPO frequency sp ectrum is simi l ar to those of so me LFQPOs from 



GRS 1915+105 during more steady states (|Ou et al. 



2010 : 



Yan et al 



2Q13). The LFQPO 



amplitude spectrum is also similar to that of a 4 Hz QPO from GRS 1915+105 during the % state 



(Zdziarski et al. 



20051) . Both of the similarities suggest that these LFQPOs from GRS 1915+105 



during different states might share the same origin. 



In summary, we have obtained detailed p cycle evolutions of LFQPO parameters. An 
elaborate comparison analysis is carried out between the results of timing and spectral analyses 
of GRS 1915+105 during p state. It shows that the LFQPO is likely originated from the corona. 
The similar spectra of variability parameters indicate that the LFQPOs of GRS 1915+105 in 
very different states seem to be produced from the same mechanism. The continuum / LFQPO 
amplitude ratio is used to probe the origin of the aperiodic X-ray variability from the accretion 
flow. 
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Fig. 1.— (a) The phase-folded PCA p class light curve in 2.0-60 keV band, (b) the 0.5-10 
Hz rms amplitude, (c) LFQPO frequency, (d) LFQPO amplitude, and (e) LFQPO coherence, as 
functions of p cycle phase. In interval I, there is no obvious LFQPO. As phase increases, the 
LFQPO frequency decreases (II and III), flattens (IV) and then increases again (V). We double 
the LFQPO frequencies in interval VI and find that they (the gray points) and the other LFQPO 
points together show a much more symmetric time evolution over the cycle. In intervals II and 
VI, the LFQPO amplitude is positively correlated with the LFQPO frequency. While in intervals 
III, IV and V, the LFQPO amplitude is negatively correlated with the LFQPO frequency. As phase 
increases, the coherence of LFQPO reveals a rapid increase (II), and then decreases quickly (III) 
followed by a continuous gentle decrease (IV, V and VI). The horizontal bars denote phase ranges. 
The vertical bars are error bars. 
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Fig. 2. — the LFQPO absolute amplitude as a function of p cycle phase. The LFQPO absolute 
amplitude increases very slightly in interval IV but dips in phase intervals 0.12-0.4 (II and III) and 
0.74-0.02 (V and VI). The horizontal bars denote phase ranges. The vertical bars are error bars. 
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Fig. 3. — (a) the LFQPO frequency spectrum in interval IV. The LFQPO frequency is positively 
correlated with the photon energy. The points are fitted with Least Squares, (b) the 0.5-64 Hz, 
3-64 Hz, 10-64 Hz continuum amplitude spectra and the LFQPO amplitude spectrum in interval 
IV. As photon energy increases, these amplitudes increase and then flatten out. (c) The ratios of 
the 0.5-64 Hz, 3-64 Hz, 10-64 Hz continuum amplitudes to the LFQPO amplitude. As photon 
energy increases, the rati os decrease rapidly and then smoo thly level off. The points can be well 



fitted with function (1) in IShaposhnikov & Titarchuk 



(120071) . The open circles, triangles, squares 



represent the 0.5-64 Hz, 3-64 Hz, 10-64 Hz continua, respectively. The filled circles represent the 
LFQPO. The horizontal bars denote energy bands. The vertical bars are error bars. 
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Fig. 4. — The ratio of the power-l aw index to th e 



The index is got from Figure 7 in Neilsen et al 



LFQ PQ frequency as a function of p cycle phase. 



(|201 lb . In interval IV, the frequency and the index 



are both relatively stable, and the ratio is ~ 0.4. As the frequency increases, the index increases, 
while the ratio decreases faster. The horizontal bars denote phase ranges. The vertical bars are 
error bars. 
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